Are stellar over-densities in dwarf galaxies the 
"smoking gun" of triaxial dark matter haloes? 
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Abstract. We use N-body simulations to study the tidal evolution of globular clusters (GCs) in 
dwarf spheroidal (dSph) galaxies. Our models adopt a cosmologically motivated scenario in which 
the dSph is approximated by a static NFW halo with a triaxial shape. For a large set of orbits 
and projection angles we examine the spatial and velocity distribution of stellar debris deposited 
during the complete disruption of stellar clusters. Our simulations show that such debris appears as 
shells, isolated clumps and elongated over-densities at low surface brightness (> 26 mag/arcsec 2 ), 
reminiscent of substructure observed in several MW dSphs. Such features arise from the triaxiality 
of the galaxy potential and do not dissolve in time. Stellar over-densities reported in several MW 
dSphs may thus be the telltale evidence of dark matter haloes being triaxial in shape. We explore 
a number of kinematical signatures that would help to validate (or falsify) this scenario. The mean 
angular momentum of the cluster debris associated with box and resonant orbits, which are absent 
in spherical potentials, is null. As a result, we show that the line-of-sight velocity distribution 
may exhibit a characteristic "double-peak" depending on the oriention of the viewing angle with 
respect to the progenitor's orbital plane. Kinematic surveys of dSphs may help to detect and identify 
substructures associated with the disruption of stellar clusters, as well as to address the shape of the 
dark matter haloes in which dSphs are embedded. 
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INTRODUCTION 

Dwarf spheroidal (dSph) galaxies are the most dark-matter (DM) dominated galaxies 
known to date (see reviews by Mateo 1998 and Gilmore et al. 2007). Dynamical con- 
straints on the DM distribution in dSphs may therefore be compared directly with the 
predictions from cosmological simulations to study the elusive nature of dark matter. 

In the last few years, however, wide- field imaging and multi-object spectroscopic 
surveys have revealed that dSph stellar components are more complex than previously 
believed. Localized oddities such as shells (Coleman et al. 2005; Olszewski et al. 2006), 
butterfly-shaped substructures (Stetson & Hesser 1998) and cold clumps (Kleyna 2003; 
Ibata et al. 2006) have been reported in several dSphs. These results conflict with the 
naive expectation that dSph DM haloes (with typical crossing times of a few hundred 
Myr) should host well-mixed stellar populations, and may thus pose a challenge to Cold 
Dark Matter (CDM). 

A second oddity in this context is the existence of globular clusters (GCs) in dSphs. 
Although GCs are considerably more common in dwarf irregular galaxies (e.g. Georgiev 
et al. 2008), at least two of the Milky Way dSphs — Fornax and Sagittarius — contain GCs 
(five and four, respectively). The formation of GCs in a cosmological context still lacks 



theoretical understanding, and it is unclear whether GCs could also have formed in other 
dSphs, but not survived to the present epoch. The present contribution elaborates on this 
scenario by studying the possible imprints that tidally disrupted clusters may leave in the 
stellar populations of dSphs. In this contribution we provide a summary of the results 
found by Penarrubia, Walker & Gilmore (2009). We refer interested readers to that paper 
for details. 

MODELS AND NUMERICAL SETUP 

Our models assume that the mass profile of a dSph can be approximated by a triaxial 
Navarro, Frenk & White (1997, NFW for short) DM halo. We choose axis-ratios c/a = 
0.5; b/a — 0.83, typically found in cosmological simulations (e.g. Jing & Suto 2002). We 
scale the free parameters of the NFW potential to those typically estimated in Milky Way 
dSphs, a peak velocity v max = 20.6 km/s and a peak radius r max = 4.1 kpc (Penarrubia 
et al. 2008a,b). 

Our GC models are N-body realizations of a King (1966) profile orbiting within a 
static NFW potential. For illustrative purposes, we choose a GC model that resembles 
the cluster Fl of the Fornax dSph, which has a mass of 3.7 x 1O 4 M , and core and tidal 
radii of R c = 10 pc and R t = 60 pc, respectively. Dynamical friction between the GC 
and the background DM particles is introduced in our computations as an external force 
term following the method outlined in Penarrubia et al. (2006). 

We follow the evolution of the cluster N-body model in the dwarf potential using 
Superb ox, a highly efficient particle-mesh gravity code. Our grid resolution is 0A6R C 
and the time step adopted is dt = t cr /25, where t CI = 15.7 Myr is the crossing time of 
our cluster model. 

RESULTS 

Triaxial potentials allow an infinite number of orbits, which are typically classified in 
four major families (i) boxes, (ii) loops (iii) resonant (iv) irregular orbits. Their phase- 
space structure show remarkable differences (e.g. Binney & Tremaine 2008). 

In this contribution we only consider clusters orbiting in the X — Z symmetry planes; 
recall that the coordinates (X, Y,Z) are oriented along the major, intermediate and minor 
axis of the potential, respectively. We emphasize that our orbital sample is a largely 
reduced subset of the whole orbital family spectrum allowed by the potential adopted in 
our study. 



Debris morphologies in triaxial potentials 

In Fig. 1 we show a few examples of tidal debris configurations that one would expect 
from a fully disrupted Fl cluster moving on orbits with different circularity (rj). 

The projected morphologies are extremely rich. By colour-coding cluster particles ac- 
cording to their projected surface brightness we can observe the formation of "folds" 
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FIGURE 1. Projected surface brightness associated with the tidal debris from the cluster Fl on different 
orbits initially confined on the symmetry plane formed by the major and minor axis of the potential. The 
initial apocentric distance is ro = 0.5 kpc and 77 is the orbital circularity. For each model we show the 
projection along the three principal axis. 



and "cusps" in regions densely populated by debris material. These over-densities 
may produce striking features. For example, folds associated with box orbits exhibit 
a "butterfly"-like shape (e.g. r\ = 0.8). Even more complex do appear the stellar sub- 
structures that arise from the disruption of clusters on resonant orbits. A good example 
is the "fish" orbit in this Figure (rj = 0.5). Note, however, the strong dependence of 
debris morphology on the line-of-sight projection. 

Our models show that stellar substructures arise naturally from the disruption of GCs 
in a static potential. Shells (folds) and isolated over-densities (cusps) arise in almost 
all orbital configurations and line-of-sight projections. These features result from the 
triaxial shape of the galaxy potential and do not dissolve in time. 



Kinematic Signatures 

The tidal disruption of a 1O 4 L globular cluster leaves stellar substructures with sur- 
face brightnesses /i^26 mag/arcsec 2 . Some of these features fall within the detection 
limits of present instrumentation. However, in areas where dwarf field stars dominate in 
number, the identification of cluster debris as stellar over-densities may prove challeng- 
ing. Kinematic surveys may provide additional constraints. 
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FIGURE 2. Projected line-of-sight velocity distributions of orbits with different apocentres (ro) mea- 
sured at the centre of the dwarf. The initial orbital circularity in these models is r\ = 0.8. 



Box, resonant and irregular orbits, all absent in spherical potentials, have null angular 
momentum. In practice this implies that, at any position angle, there is the same amount 
of debris material moving away from us, as towards us. This results in a line-of-sight 
velocity distribution that is always symmetric and centred at (v/ 5 ) = 0, but which may 
show distinctive "double -peaks" depending on the oriention of the viewing angle with 
respect to the progenitor's orbital plane, as shown in Fig. 2. Such velocity distributions 
would not be expected in well-mixed stellar populations. Kinematic surveys of dSphs are 
thus useful tools to identify substructures associated to the disruption of stellar clusters, 
as well as to address the shape of DM haloes. 
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